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Inspired from an experimentally examined organic conductor, a novel mechanism for negative 
magnetoresistance is proposed for repulsively interacting electrons on a lattice whose band dispersion 
contains a flat portion (a flat bottom below a dispersive part here) . When the Fermi level lies in the 
flat part, the electron correlation should cause ferromagnetic spin fluctuations to develop with an 
enhanced susceptibility. A relatively small magnetic field will then shift the majority-spin Fermi level 
to the dispersive part, resulting in a negative magnetoresistance. We have actually confirmed the 
idea by calculating the conductivity in magnetic flelds, with the fluctuation exchange approximation, 
for the repulsive Hubbard model on a square lattice having a large second nearest-neighbor hopping. 
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Negative magnetoresistance has provided fascination 
in diverse classes of systems, such as impurity band in 
semiconductors or Mn oxides. The problem is the in- 
terplay of the transport and the spin structure, and how 
the spins and carriers respond to magnetic fields can vary 
from system to system. In the physics of correlated elec- 
tron systems, spin is a key ingredient. In fact, itiner- 
ant (metallic) ferromagnctism has beeri_a central problem 
from the days of Kanamoritl, HubbardH and GutzwillerO. 
It is becoming increasingly clear that the criterion for 
the ferromagnetism becomes stringent as one incorpo- 
rates correlation effects, but that a large density of states 
at the Fermi level does favor ferromagnetism. 

In this paper, we propose a novel mechanism for neg- 
ative magnetoresistance, which is realized as a combined 
effect of the electron correlation (proximity to ferromag- 
netism) and the band structure (coexistence of flat and 
dispersive parts). The idea is the following: We con- 
sider repulsively interacting electrons on a band, whose 
one-electron dispersion has a flat part in an otherwise 
dispersive band. For low electron densities with the 
Fermi level Ep located near the flat part of the band, the 
system should have a strong tendency toward ferromag- 
netism. A large density of states D{Ep) at the bottom 
of the band is indeed a situation originally considered by 
Kanamoru and received a renewed ^interest recently for 
two and three dimejiAional systemfflU as well as in one di- 
mensional systemaJfl. When a magnetic field is applied, 
a weak field is then enough to drive the system into a 
significant magnetization (spin polarization), since the 
magnetic susceptibility is enhanced due to ferromagnetic 
fluctuations. The Ep for the majority spins will then 
be shifted into the dispersive (i.e., lighter-mass) part of 
the band, and we expect that the system becomes more 
conductive for relatively weak magnetic fields. When the 
polarization is small enough there is a possibility that 
the effect of the increase in vp of the majority spin may 
be compensated by the decrease in vp of the minority 
spin, but such a compensation will not occur when the 



polarization is sufficiently large due to the enhanced sus- 
ceptibility, since not only vp but also the number of the 
majority spin increase in the presence of the magnetic 
field. Further, there will be fewer electron-electron scat- 
terings for a polarized state, at least for the Hubbard 
model where only up and down spin electrons interact. 
These will make the system more conductive. This idea 
has been inspired from a recent experimental result on a 
certain class organic conductors, called r-type conduc- 
tors, for which Murata et aB have observed negative 
mangetoresistance, see below. The band structure cal- 
culation indicates that the dispersion is indeed flat at 
the bottom along and ky directions. 

In the following, we confirm this idea for the single- 
band Hubbard model, a simplest model for repulsive 
electron correlation, on a simplest flat-bottomed tight- 
binding model, the square lattice having a large second 
nearest neighbor hopping, t', along with t between near- 
est neighbors. The one-electron dispersion becomes 



2t (cos kx + cos ky) + At' cos k^ cos ky 



(1) 



and the dispersion is flat along k^ and ky (right panel 
of Fig. |l]) for t' ~ 0.5|t| with the van Hove singularity 
coming down to the bottom. 

For t' ~ 0.5 (we take |t| = 1 hereafter), Hlubina has 
shown, with the T-matrix approximation for the Hub- 
bard modeljnjthat the ground state is fully spin-polarized 
for n ~ 0.4.QO Since the T-matrix approximation is valid 
only in the low enough densities and for small U , here 
we start with obtaining the diagram for large U with the 
exact diagonalization of a finite system (4x4 sites with 8 
electrons with the band quarter filled). 

Detection of ferromagnetic states from total in fi- 
nite systems requires_some care as we have revealed in 
previous publicationscl. Namely, ferromagnetic states in 
itinerant electron systems can accompany a spiral spin 
state in finite systems. The spiral state is a spin singlet 
state having the spin correlation wave length as large as 
the linear dimension of the system, which can lie in en- 
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ergy below the ferromagnetic state. We found that the 
ground state of the 16-site system has always = 
in periodic boundary conditions in both x and y, but if 
we look at the spin structure iS'(q) (Fourier transform 
of {Si ■ Sj)) the state turns out to be indeed spiral (i.e., 
>S'(q) peaked at (0, ztir/L) and (±7r/i, 0)) for large U and 
f ~ 0.5. 

Although we can identify the spiral state as ferroma^ 
netic in the thermodynamic limit, it has been shownE2l 
that a faster approach to the thermodynamic limit is at- 
tained if we adopt an appropriate boundary condition 
(periodic in one direction and antiperiodic in another) 
to selectively push the spiral state above the ferromag- 
netic state in energy. If we look at the phase diagram 
thus obtained against t' and U in Fig. |l|, the fully polar- 
ized ferromagnetic region is seen to exist for t' ~ 0.5 and 
U >6. 

Keeping this phase diagram in mind, we now discuss 
the conductivity in magnetic fields B. The conductiv- 
ity of the interacting system is calculated here by the 
fluctuation exchange approximation (FLEX). The FLEX, 
introduced by Bickers et aLtil, treats spin and charge 
fluctuations by starting from a set of skeleton diagrams 
for the Luttinger-W^d functional, based on the idea of 
Baym and KadanofH. Then a (fc-dependent) self energy 
is computed from RPA-type bubble and ladder diagrams 
self-consistently. 

The dc conductivity is given in the Kubo formula as 



7^.= lime^^ 
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ImXkk'crcr' (W + iS) 



where -ftTkk' o-o-' (w + iS) is the Fourier component of the re- 
tarded two-particle Green's function and v^^ — 
is the unperturbed velocity. rWe set h = 1, fcs = 1 here- 
after. If we follow Eliashbergt3 for the analytic continua- 
tion of K(uj + iS) from K{iLUn), where Un is the Matsub- 
ara frequency, the conductivity per spin reads, for long 
enough life time of the quasi-particle. 
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where v^^^ is the dressed velocity, J^^ the current, 7k the 
damping constant of the quasiparticle, and / the Fermi 
distribution function. 

Note that the conductivity of a clean system of inter- 
acting electrons can be finite even though the interaction 
is an internal force, since the momentum is dissipated 
through UmklapjjHprocesses. It has been shown by Ya- 
mada and Yosidall3 that the conductivity diverges, as it 
should, when the Umklapp processes are turned off only 
if we consider the vertex correction for the current Jk 
appropriately, so the correction should be included for a 
consistent treatment. Within the FLEX, which is a con- 
serving approximation—there are three types of diagrams 
for irreducible vertextJ, namely two Aslamasov-Larkin 
(AL) type diagrams and one Maki-Thompson (MT) type 



one. Kontani et alii have shown that, when antiferro- 
magnetic fiuctuations are dominant, the AL contribution 
can be neglected. We can extend this argument to show 
that the AL term is also negligible when ferromagnetic 
fluctuatioas. are dominant. Hence we consider only the 
MT term.E3 

Let us now discuss the conductivity in magnetic 
flelds. We have then to add the Zeeman term, 
/i ^]^.^ sgn(cr)C]{^g.CkCT, to the Hamiltonian, where h = 
g^isB is the Zeeman energy with <? ~ 2. Green's function 
and other quantities then become ti-dependent. To con- 
centrate on the effect of the Zeeman splitting we assume 
here that the direction of the magnetic field is parallel to 
the current, so that we do not have to take account of 
the effect of the field on orbital motions. 

The conductivity is obtained from the Bethe-Salpeter 
equation, which is a simple extension of the equatiops 
derived by Kontani et al for the spin-independent casell3 
to the present spin-dependent case. We end up with the 
diagonal conductivity. 
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where N is the number of sites, T the temperature, 
Gko-(ijj) the dressed Green's function, and the velocity 
Vkxa = (9/9fcx)[£k + Rei;kcr(t^ = 0)] with Skcr(a^) being 
the self energy. 

The kernel, V\^aa' (^), which contains the effect of fluc- 
tuation exchanges, caiL_be obtained by an analytic con- 
tinuation of Vkcrcr' (iti'ri)t-J, 
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where k = (k, iujn) in the last line. 

Let us now present the results. We take the case of 
t' = 0.5, U = 2 with the band filling n = 0.4, which falls 
upon the paramagnetic region close to the ferromagnetic 
boundary in the phase diagram, Fig.|l|. We first check 
that we do have strong ferromagnetic fluctuations. Fig- 
ure ^ shows the wave number dependence of the spin 
susceptibility, Xk^^ = 2xl^(w = 0)/(l - UxU^ = 0)) for 
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T = 0.1. We can see that there is indeed a peak around 

r(k = o). 

Let us next focus on the static, uniform magnetic sus- 
ceptibihty. Since we are sitting close to the ferromag- 
netic boundary, the susceptibihty should be finite but 
enhanced. Here we compute the quantity in two ways: 
one is to calculate the derivative x = d{na- — n-a)/dh 
with a finite difference 6h = 0.005, and the other is Xn^^- 
If we look at their temperature dependence in Fig. for 
t' = 0.5, n = 0.4 with various 0.5 < i7 < 2.0, they both 
sharply increase for T — > 0. There is a deviation between 
X and Xo^^ for larger U. The deviation itself indicates 
that the irreducible four-point vertex T — 
cannot be approximated by U. Since we are dealing with 
a case where the ferromagnetic phase appears for larger 
U{> 6 ~ 7), X underestimates the effect of external mag- 
netic field for large U. 

Now, we come to the key result for the magnetoresis- 
tance, given in Fig. | for t' = 0.5, U = 2, n ^ 0.4. The 
figure compares the resistivity against T in the absence 
and in the presence of the magnetic field. We take the 
system size A'' = 64^, with 512 Matsubara frequencies, 
which is checked to be sufficient for the temperature re- 
gion studied here. We can see that we do have a negative 
magnetoresistance. The change in the resistance is of the 
order of 10% for the Zeeman energy h = 0.05. If we were 
not sitting close to the ferromagnetism, much larger fields 
would be required. The vertex correction does not alter 
the result significantly, which implies that Boltzmann's 
transport picture (with no vertex correction) already has 
the negative magnetoresistance. The negative magne- 
toresistance should become more prominent at lower tem- 
peratures where the ferromagnetic fluctuation increases. 
This is not too noticeable in the result, which should be 
an artifact: as mentioned above, the spin polarization 
is underestimated in the FLEX at low temperatures for 
[/ = 2. 

As tooiched upon at the beginning, Papavassil- 
iouet aZJla has found experimentally that organic salts 
DaAiAy, based on D (= P-S, S'-DMEDT-TTF or EDO- 
S',S'-DMEDT-TTF) in combination with linear anions A 
(=AuBr2, I3, or IBr2) with the fractional value of y con- 
trolling the density of carriers are two-dimensional met- 
als in the r crystal form. The band structure of a single 
layer of the r-phase, calculated with tha^xtended Hiickel 
method, contains a flat-bottomed bando. 

In terms of the tight-binding model, we can regard 
that the in-plane molecular configuration is such that 
the next-nearest intermolecular hopping, t2 — 0.02 eV, 
appears in every other plaquetts in a checker-boatd man- 
ner on top of the nearest one ti ~ 0.2 eV (Fig. ||)E3. The 
checker-board makes the Brillouin zone folded, where 
the dispersion of the upper band in which Ep resides 
is like the square root of Eqn.(l) with t' = 0.5 when the 
splitting due to <2 is small. This way we can have an 
anisotropically flat dispersion along , ky with a large 
D{Ep) at the bottom of the band, so a ferromagnetic 
component in the spin fluctuation exists for this disper- 



sion, as we have checked with FLEX. If we naively plug in 
1/8 (~ 0.05 eV) of the width of the upper band calculated 
in ref. |^ for i, the h = 0.05 for [/ = 2 in Fig. | corre- 
sponds to S ~ 30 T. Since the divergence of D{Ef) is 
weaker than that for Eqn.(l) because of the square root, 
large values of U should be necessary to obtain apprecia- 
ble negative magnetoresistance. In organic materials U 
may be indeed large, but we have then to check whether 
the band mixing across the gap for large U can smear 
out the ferromagnetism. Our preliminary result for the 
two-band model with the multi-band FLEX shows that 
the ferromagnetic component in the spin fluctuation does 
remain for J7 ~ 4, although it becomes weaker. Thus 
quantitative estimates of the magnetoresitance for the 
material should include these effects. 

Murata et am have in fact observed a negative magne- 
toresistance for B ~ sevecal tesla with some hysteretic 
behaviors in this materialEll. Theoretically, however, we 
can argue that dominant ferromagnetic fluctuations are 
not a necessary condition for the negative magnetore- 
sistance conceived here. In the present scheme, the the 
system becomes more conductive in moderate B due to 
a combination of a correlation-enhanced spin suscepti- 
bility and the flat-bottomed dispersion. So all we have 
to have about the magnetism is a large enough suscepti- 
bility. Thus the negative magnetoresistance observed by 
Murata et aEI is understandable provided the ferromag- 
netic component is present, if not dominant, in the spin 
fluctuation for the present mechanism to be relevant. 

Although we have exemplifled our idea so far in two di- 
mensions, we believe that the present mechanism should 
be general, and can be found in other models with strong 
ferromagnetic spin fluctuations such as face centered cu- 
bic lattice, whose band structure has flat and dispersive 
parts as well. 
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A preliminary result (K. Murata, private communications) 
indicates both an enhanced susceptibility from an M — H 
curve and antiferromagnetic fluctuations from NMR. The 
present model for r-conductor should in fact be strongly 
antiferromagnetic if t2 were absent, since the system re- 
duces to a nearly half-filled square lattice. On the other 
hand ferromagnetic fluctuation is dominant if the lower 
band is neglected, so antiferromagnetic and ferromagnetic 
fluctuations should coexist for finite values of t2, which will 
imply a coexistence of antiferromagnetism and the negative 
magnetoresistance. 
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FIG. 2. The wave number dependence of Xk^^ for 
t' = 0.5, C/ = 2,n = 0.4,T = 0.1. 
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FIG. 3. The temperature dependence of the static mag- 
netic susceptibility x^'^^'^ (dashed line) and Xk=o (f^U line) for 
U = 0.5, 1.0, 2.0 with t' = 0.5, n = 0.4. 
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FIG. 1. The phase diagram of the t-t' Hubbard model, de- 
termined by the exact diagonalization of 8 electrons in 4 x 
4 system. The right panel is the one-electron dispersion for 
t = -l,t' = 0.5. 
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FIG. 4. The temperature dependence of the resistivity of 
the t-t' Hubbard model with (bottom) or without (top) vertex 
corrections for t' = 0.5, ?7 = 2,n = 0.4 in zero magnetic field 
{h — 0, full line) and in a magnetic field {h — 0.05, dashed 
line) . 
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FIG. 5. The in-plane molecular configuration in 

T-D2A1AJ;. The solid lines denote the nearest neighbor hop- 
ping while dashed lines the second nearest neighbor hopping 
between face-to-face molecules. 
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